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We consider an Einstein-Scalar-Gauss-Bonnet gravitational theory, and argue that at early times 
the Ricci scalar can be safely ignored. We then demonstrate that the pure scalar-Gauss-Bonnet 
theory, with a quadratic coupling function, naturally supports inflationary - de Sitter - solutions. 
During inflation, the scalar held decays exponentially and its effective potential remains always 
bounded. The theory contains also solutions where these de Sitter phases possess a natural exit 
mechanism and are replaced by linearly expanding - Milne - phases. 
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I. INTRODUCTION 

Einstein’s theory of gravity, based on General Rel¬ 
ativity, is known to possess many successes but also 
many open questions. In the context of Cosmology, 
the “old problems” of flatness, horizon, monopoles 
and density perturbations were solved by the intro¬ 
duction of inflation, a very fast expanding, early phase 
in the history of our universe [1H3|. New questions 
have however emerged, such as the nature of dark 
matter and dark energy or the coincidence problem, 
which are still persisting despite all the efforts. 

Searching for alternative and maybe more funda¬ 
mental theories of Gravity, that would allow for new 
directions of thinking and for novel solutions, has been 
a path followed by many scientists. Looking at higher 
dimensions for inspiration, the Lovelock theory Q 
emerges as a natural higher-dimensional (d > 4) gen¬ 
eralization: its action is a homogeneous polynomial of 
degree N of Riemann curvature, with Einstein’s term 
R arising for N = 1, the quadratic Gauss-Bonnet term 
R 2 GB = R tlvp<J W lv P a - 4 R pv R» v + R 2 for N = 2, and 
so on. The most remarkable and distinguishing basic 
feature of the theory is that the ensuing field equations 
are of second order as in Einstein’s gravity, a require¬ 
ment necessary to ward off undesirable features like 
ghosts. 

Our primary candidate for a generalised gravita¬ 
tional theory in four dimensions is the onegiven by the 
heterotic superstring effective action [a-Q, where the 
Gauss-Bonnet (GB) term appears as an a'- 0rc j er C or- 
rection. In addition, this term is the highest, non-zero 


one in the Lovelock polynomial in d = 4 - in d dimen¬ 
sions, only the terms with N < d/2 are non-trivial Q. 
But being a topological invariant in 4-dimensions, the 
GB term, in the context of any theory, should always 
be coupled to a scalar field to get a non-trivial con¬ 
tribution to the field equations. The presence of the 
GB term in the string effective theory, or in more gen¬ 
eralised string-inspired ones, has been shown to lead 
to novel solutions such as singularity-free cosmologi¬ 
cal solutions [gj, H , hairy black holes mm or even 
traversable wormholes [12|. 

In all cases, the GB term was shown to be the domi¬ 
nant agent that allowed for the existence of these solu¬ 
tions. Motivated by this, here we consider an Einstein- 
scalar-GB theory with a polynomial coupling function 
between the scalar field and the GB term. We focus 
our analysis on the early-time evolution of the universe 
and argue that the Ricci scalar is again sub-dominant 
to the GB term, and thus it can be ignored. We then 
demonstrate that analytical, elegant, inflationary so¬ 
lutions with attractive characteristics - compared to 
the existing models - are naturally supported by the 
field equations. 


II. THE EINSTEIN-SCALAR-GB THEORY 


We consider the following gravitational theory of a 
scalar field <j> 
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coupled non-minimally to gravity via the Gauss- 
Bonnet term Rq B by a general coupling function /(</>). 
The variation of the action © with respect to the met¬ 
ric tensor and scalar field leads to the field equations 
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V^g 

where /' = df /d(j), the constant k 2 = 8ttG has been 
set to unity, and P^av /3 is defined as 

PfiauP RfiavP P Zg^pR^a 4” 2 g a [i/Rp](i P Rgp[i>g /3]a • 

( 4 ) 

We assume that the line-element has the 
Friedmann-Lemaitre-Robertson-Walker form 


ds 2 = —dt 2 + a 2 [t] 


dr 2 


1 — kr 2 


P r 2 (d0 2 + sin 2 6dip 2 ) 


( 5 ) 

with k = 0, ±1 denoting the curvature of the 3- 
dimensional space. Then, the field equations take the 
explicit form 
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0p3H<fi-3f'(H 2 + — )(H 2 + H) = 0, 


( 8 ) 


where H = a/a is the Hubble parameter and the dot 
denotes a derivative with respect to time. In this 
work, we will consider a polynomial coupling func¬ 
tion, i.e. f(cj>) = X (f> n , where A is a constant and n an 
integer. The case with n = 0 is equivalent to ignoring 
the GB term; the cases with n = 1 and n = 2 will be 
studied in Sections III and IV, respectively. 


III. THE ROLE OF THE RICCI SCALAR 

In this section, we will formulate an analytical ar¬ 
gument to demonstrate that, in the strong-gravity 
regime, the presence of the Ricci scalar significantly 
complicates the derivation of a cosmological solution 
without adding any important effect to its dynamics. 
Here, we will consider the case of n = 1, for which the 
scalar equation JS]) can be integrated once to yield the 
relation 


D 


Aa (3 k + a 2 ) 


with D an integration constant. The second term in 
the above relation is clearly the one associated to the 
GB term, and we will focus on this, setting D = 0. 

Combining Eqs. ©-© and using that / = Xcp, we 
obtain the constraint 

4 (H 2 P 4) + 2 (H 2 PH) (IP A H<j>) 
a z 

P3X 2 (H 2 + H)(H 2 + 4) 2 = 0 ■ (10) 

If we substitute <p via Eq. m, the above equation has 
no dependence on the scalar field. After some algebra, 
it can be integrated once to give 

4 2 f\3ka 2 5 a 4 2k 3 A . . 

3 “=- A (— + - + F^J +C ‘- (11) 

with ci an integration constant again. For k ^ 0, the 
above cannot be easily solved for d. However, for a 
flat universe (k = 0), we easily find that 
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The solution of the above equation can be written in 
the non-explicit form 


i(t)F 


115 3 a 4 (t) 
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where F(a, b , c; x) is the hypergeometric function. The 
above describes an increasing scale factor with time 
that is however bounded from above due to the conver¬ 
gence criterion a(t) < a max = (ci/3) 1 / 4 of F. In the 
limit a —> 0, i.e. in the strong-gravity limit, the rela¬ 
tion between the scale factor and the time coordinate 
becomes linear with a Big-Bang singularity emerging 
at a finite value of time. 

Let us now focus on the strong-gravity regime from 
the beginning and assume that, in that limit, the Ricci 
scalar R will be subdominant to the GB term. In this 
case, the scalar equation ® remains the same, while 
the two gravitational equations will be modified by 
loosing all unity terms appearing inside brackets on 
their left-hand-sides. Then, Eq. © can be solved for 
<j), even for a general coupling function f(<f>), to give 


, = 6 f H (H 2 + -). 


(14) 


Following a similar analysis as above, i.e. combining 
Eqs. ©-©, employing Eq. (fill) and the fact that 
/ = X(f>, we now arrive at the simpler constraint 

(5 H 2 + ^(H 2 +H) = 0 . 
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The general solution of the above differential equation 
is: a(t) = At + B , describing a linearly-expanding uni¬ 
verse at early times. This is in full agreement with the 
result obtained above in the context of the full anal¬ 
ysis and in the early-time limit. By plotting the two 
solutions [13, one may easily see that these perfectly 
match at early times, while a deviation begins to ap¬ 
pear, for a given A, the sooner the larger the value of 
the integration constant ci is. A similar conclusion 
may be reached for the behaviour of the scalar field: 
using either Eq. © or (fill) and taking the limit of 
early times, we obtain (f> ~ (At + B )~ 2 with a diver¬ 
gence appearing at the initial singularity. 


IV. INFLATIONARY EXPANSION 

In this section, we will proceed to study the case of a 
quadratic coupling function, / = A (ft 2 . Unfortunately, 
in this case, the set of equations is significantly more 
complicated and the formulation of a similar analytic 
argument regarding the role of the Ricci scalar is im¬ 
possible. However, since the coupling function deter¬ 
mines the weight of the GB term in the theory, we ex¬ 
pect that its exact form merely determines the point in 
time where the GB term begins to dominate over the 
Ricci scalar in the early-time, strong-gravity regime. 
Thus, here we will ignore again the Ricci scalar and 
focus on an appropriately chosen early-time period. 
Equation (fTdl) still holds, however, the different form 
of the coupling function modifies the constraint (1151) 
as follows 

(5 H 2 + ^)(H 2 +H) + 24A H 2 [h 2 + A) 2 = 0 . (16) 

The above differential equation involves only the scale 
factor a(t), however it is only for k = 0 that it can 
be integrated twice to yield an explicit solution - for 
k ^ 0, a transcendental equation arises after the first 
integration [l3[. Thus, for a flat space, we have 

H + H 2 ( 1-^=0, (17) 

H dS 

where we have defined H 2 S = — 5/24A. The above im¬ 
plies that a de Sitter solution, with H a constant and 
thus H = 0, exists when A < 0 describing an inflation¬ 
ary phase. The corresponding accelerating solution 
with H 2 s = —5/24A is a repeller, while the second so¬ 
lution H = 0, describing a Minkowski spacetime, is an 
attractor. Therefore, in order to have a natural exit 
from the inflationary phase we have to consider an 
initial condition (Hi) such that Hi < Hds■ Note that 
if we consider Hi > Hds, the Universe will be eter¬ 
nally inflating (H > 0). Integrating Eq. ( 1171 ) once, we 



FIG. 1: The scale factor as a function of time for Ci = 0 
and |A| = 0.1, 0.2, 0.5,1, 2,10 (from top to bottom). 


obtain 


H = 


H dS 



2 


(18) 


where Ci is an arbitrary integration constant. De¬ 
pending on the values of the two parameters A and 
Ci, the above equation results to different types of cos¬ 
mological solutions. The existence of de Sitter implies 
A < 0 and the natural exit implies Ci > 0 (H < Hds)- 
Here, we will present only these two cases that support 
a non-eternal inflationary phase of expansion for the 
universe - for a more comprehensive analysis, please 
refer to [bit . 


A. The case with Ci = 0 and A < 0 


If Ci = 0, then integrating once more Eq. m , we 
find 


a(t) =a 0 exp 2 ^a[ ' ( 19 ) 

The profile of the scale factor is shown in Fig. 1. 
We observe that the smaller the value of the coupling 
constant |A| is, the faster a(t) evolves with time. The 
above solution may be considered as an alternative 
to the usual inflationary solutions following from the 
Einstein-Hilbert action which require an appropriate 
potential for the scalar field here, it is the GB term 
itself that provides the necessary potential. 

The solution for the scalar field can be found via 
Eq. (Il4l) : for k = 0 and f = 2A </>, it can readily be 
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FIG. 2: Cosmological solutions with C\ > 0, A < 0, and 
for v = 0.1, 1,10,100. 


integrated to give the solution 

m = exp ■ (20) 

The above describes a regular, exponentially decaying 
scalar field. One may easily verify that the combina¬ 
tion of solutions (EH) and (l20l) satisfies the complete 
set of field equations. 

The scalar field starts from an arbitrary value (f) o 
and quickly reduces to zero - again, the smaller the 
value of | A| is, the faster it decays. As a result, the cou¬ 
pling function f((j>) = A cj) 2 remains always bounded. 
The effective potential of the scalar field receives con¬ 
tributions from both the GB term and the coupling 
function and has the explicit form 

I 2^ (h 2 

/«)/&, = ^. (2i) 

Therefore, V e ff also remains bounded as the field 
evolves; at the same time, it may become arbitrar¬ 
ily large at early times by appropriately choosing the 
value of the coupling constant A. This is in sharp con¬ 
trast to the behavior seen in more traditional inflation¬ 
ary models, such as the chaotic (H, where the effective 
potential blows up for the super-Planckian initial val¬ 
ues of the field unless a fine-tuning is imposed on the 
parameters of the model - as we will see in Sec. IV, 
such an unnatural initial value is not needed in our 
case for the necessary number of e-foldings to be ob¬ 
tained. It is also a significant improvement compared 
to models for inflation where the scalar potential 
remains bounded but its value is too small to justify 
its dominance over any other distribution of matter in 
the universe and thus allow inflation to be realised. 


B. The case with C i > 0 and A < 0 


We will now demonstrate that the pure de Sitter 
solution of the previous case is the early-time limit 
of a more interesting class of cosmological solutions 
that follow when A < 0 and C\ positive. In this case, 
by using the change of variable a = v taniu, where 
v 2 = 12|A|/Ci, Eq. (1TH1) can be again integrated to 
yield 


r-z - - ( \/a 2 + v 2 — v \ / 5 . . 

+ ^‘ n (—;—j = V2cr ( o) ' 

( 22 ) 

Plotting the above (see Fig. 2), the scale factor is 
found to be an increasing function of time. In the limit 
a —> 0, i.e. at early times, Eq. (l22ll reduces to the pure 
de Sitter solution (fl9l) : on the other hand, expanding 
P^l) for a 2 3> u 2 , we obtain a linear dependence of the 
scale factor on time, a{t) ~ ^/5/2Ci t. Thus, in this 
case the universe naturally interpolates between a de 
Sitter phase at early times and a Milne phase at later 
times. 

In order to find the solution for (j>, we observe that, 
for /(</>) = A cj) 2 , it holds that f" = Then, using 

Eq. (fl4l) in Eq. (fTfil) and integrating once, we obtain 
the relation </> 2 a 5 = const. = Co- Then, the following 
implicit expression may be derived for the scalar field 


,- 0 ( ¥) .- s£ ^ 


(23) 


For early times, (f> is again decreasing exponentially 
while, at later times, it reduces to a constant. Its ef¬ 
fective potential is also decreasing, exponentially at 
early times and as 0(t~ 6 ) at later times, remaining 
again bounded. However, we expect that the late-time 
behaviour of both the scale factor and the scalar field 
will start to deviate from the above expressions due to 
the increasing role of the Ricci scalar as the universe 
expands and to the contribution of radiation/matter. 
Only a numerical analysis of the full model could re¬ 
veal the duration of the Milne phase and its impli¬ 
cations for late cosmology as well as the fate of the 
scalar field. 

The early de Sitter phase should have a long-enough 
duration in order to resolve the problems of the stan¬ 
dard cosmology. Using the form (fTSl) . we demand the 
initial condition to satisfy Hi ~ Hds{ 1 — e~ 2N ) which 
implies C\ ~ \\\e~ 2N , where N is the number of e- 
foldings (N > 60). According to this, the smaller the 
value of Ci, the longer the acceleration; this was ex¬ 
pected since for Ci = 0, the acceleration is eternal. 
Therefore the number of e-foldings during inflation 
depends on the initial condition Hi. 
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FIG. 3: Evolution of the slow-roll parameter e\ as function 
of N = log(a) for different initial conditions Ci = \\\e~ 2N 
which gives an e-folding of (40,60,80), respectively. 


We may finally define the slow-roll parameters ei = 
—H/H 2 and ej+i = —ii/Ha as follows 

e 2n+l = d = 1 - H 2 /H 2 dS , (24) 

e 2 „ = e 2 = 2H 2 /Hj s . (25) 

During inflation, the odd slow-roll parameters ap¬ 
proach zero, C 2 n-i — 0 , while the even parameters 
are of order unity, 62 n — 2. The slow-roll approxima¬ 
tion is thus violated. A very similar behavior exists in 
the so called ultra-slow roll inflation 0 (also known 
as fast-roll inflation 0 )- 

At later times, the slow-roll parameters move away 
from the above values to mark the end of inflation. 
In Fig. [3] we depict the evolution of ei, for different 
initial conditions: inflation is realized at early times, 
when ei ~ 0 , and a natural exit takes place at later 
times as ei moves to values closer to unity. 

Let us finally note that the cases of a general poly¬ 
nomial coupling function, i.e. with /(</>) = A</>" and 
n an arbitrary integer, as well as that of an exponen¬ 
tial coupling function, i.e. with f(cf>) = Ae K ^ and k 
an arbitrary number, were also studied fl3| . In the 
first case, it was demonstrated that de Sitter-type 
solutions, either as exact or limiting solutions, arise 
only for the case of n = 2. In the latter case, it was 
shown that a scalar-GB theory with an exponential 
coupling function admits, at early times, only slowly 
expanding solutions of the form a(t) ~ (At + B) 1//5 . 
Therefore, the particular model apparently singles out 
the quadratic coupling function as the unique choice 
that allows for the emergence of inflationary solutions 
or solutions with additional attractive characteristics 
from the cosmological point of view 0 ] . 


V. CONCLUSIONS 

We have considered an Einstein-Scalar-GB gravita¬ 
tional theory, and developed an analytical argument 
to show that the Ricci scalar can be safely ignored 
when we focus on early times. Note that in the high 
energy regime of the early Universe, higher-order cur¬ 
vature terms would naturally be the most pertinent. 
Here, we have considered only those terms inspired 
by the Lovelock expansion, that lead to second-order 
differential equations, a necessary condition to ward 
off ghost fields from the theory. As the higher Love¬ 
lock polynomials, with N > 3, are trivial in four di¬ 
mensions, even if coupled to a scalar field (since the 
Lovelock Lagrangian vanishes, C N = 0, for N > 3), 
it is only the GB term that could have a non-trivial 
contribution, through a scalar field coupling, on the 
4-dimensional dynamics. 

We have demonstrated that a pure scalar-GB the¬ 
ory naturally supports de Sitter solutions with a nat¬ 
ural exit mechanism to a linearly expanding - Milne 
- phase. The quadratic coupling function, although a 
special choice, is in fact the only one that allows for 
inflationary, de Sitter-type solutions to emerge [l3j; 
in addition, it shares several features with the one 
emerging in the context of the heterotic superstring 
effective action and was shown in the past to support 
similar type of singularity-free cosmological solutions 
[fl]. Our simplifying approach leads to analytical, el¬ 
egant inflationary solutions with a bounded potential 
and free of fine-tunings. Important questions that still 
need to be answered are how, and at which point in 
time, the Ricci scalar and the presence of additional 
matter affect the dynamics of the model - for this, a 
more involved numerical analysis is clearly necessary. 
Also, one should investigate whether this type of in¬ 
flation gives a perturbation spectrum in accordance 
to observations. Note, however, that all the afore¬ 
mentioned questions will be relevant not only to our 
model but to any string-inspired or generalised gravi¬ 
tational theory that contains the Ricci scalar and GB 
term. We plan to report on these questions in a future 
work. 
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